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In the Cadassian universe, one can explain the acceleration of the universe without introducing 
dark energy component. However, the dynamical equations of this model can not be directly ob- 
tained from the action principle. Recently, works on the relation between thermodynamics and 
gravity indicates that gravity force may not be the fundamental force. In this paper, we study the 
thermodynamics of the Cardassian universe, and regard it as the origin of this cosmological model. 
We find that the corresponding entropy obeys ordinary area law when the area of the trapping hori- 
zon is small, while it becomes a constant when area is going to be large in the original and modified 
polytropic Cardassian model, and it has a maximum value in the exponential one. It seems that the 
Cardassian universe only contains finite information according to the holographic principle, which 
states that all the information in the bulk should be encoded in the boundary of the bulk. 

PACS numbers: 



I. INTRODUCTION 

Recently, the relation between the laws of thermody- 
namics and that of gravity has aroused great interest. In 
Hawking's work [H, it shows that black holes have ther- 
mal radiation with a temperature determined by its sur- 
face gravity due to the quantum effect, which implies that 
these two branches of physics may have underly relations, 
namely, they are unified. The authors in ref. [2| proposed 
a unified first law of black hole dynamics and relativistic 
thermodynamics in spherically symmetric general rela- 
tivity, see also By using the unified law in the 
framework of Friedmann-Robertson- Walker (FRW) uni- 
verse, the authors in 6] find the relation between the 
Friedmann equation and the first law of thermodynamics 
on the "inner" trapping horizon. They also find that the 
Clausius relation holds for the Gauss-Bonnet and Love- 
lock theory by treating the higher derivative terms as 
an effective energy-momentum tensor, but the Clausius 
relation never holds for the scalar-tensor theory, which 
implies that this is a system with non-equilibrium ther- 
modynamics. These results indicate that the gravity may 
be not one of the fundamental forces and it is only an- 
other view of thermodynamics. In this spirit, we can 
regard the thermodynamics law as the fundamental law 
and it will give the corresponding gravity theory. Thus, 
one can study the thermodynamics of some cosmological 
models to pursue their origins. For related discussions, 
see also |7H13j. 

Up to now, there are many kinds of dark energy mod- 
els and modified gravity theories proposed to explain the 
current accelerating expansion of the universe, which has 
been confirmed by the observations like Type la super- 
novae (SNe la), CMB and SDSS et al. The dark en- 
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ergy models assume the existence of an energy compo- 
nent with negative pressure in the universe, and it dom- 
inates and accelerates the universe at late times. The 
cosmological constant seems the best candidate of dark 
energy, but it suffers the fine tuning problem and coin- 
cidence problem, and it may even have the age problem 
p^ . To alleviate these problems, many dynamic dark 
energy models were proposed. However, people still do 
not know what is dark energy. 

Since the Einstein general gravity theory has not been 
checked in a very large scale, then one does not know 
whether this gravity theory is suitable or not for studying 
the observational data like SNe la, and maybe the accel- 
erating expansion of universe is due to the gravity theory 
that differs from the general gravity. Thus, many modi- 
fied gravity theories like /(i?), DGP et al. are proposed 
to explain the accelerating phenomenology. The Cardas- 
sian model is a kind of model in which the Friedmann 
equation is modified by the introduction of an additional 
nonlinear term of energy density and in this model one 
does not need dark energy component also. However, one 
can not directly find its origin from the first principle. 

As we mentioned before, the gravity may be not a fun- 
damental theory, and one can at least study the thermo- 
dynamics of a cosmological model to pursue its origin. 
So, in this paper, we will study the thermodynamics of 
the Cardassian universe, and regarded it as its origin. 
We find that the corresponding entropy obeys ordinary 
area law S = A /{AG) on the trapping horizon when A is 
small, where A is the the area of the trapping horizon, 
and it becomes a constant when A is going to be large in 
the original (OC) [15] and modified polytropic Cardas- 
sian model (MFC) [16], while it has a maximum value 
in the exponential model (EC) [U [ill- It seems that 
the Cardassian universe only contains finite information 
according to the holographic principle, which states that 
all the information in the bulk should be encoded in the 
boundary of the bulk. 

This paper is organized as follows: In Section [TTl we 
will give a briefiy review on the unified first law and its 
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application. In Section IIIIl we will study the thermo- 
dynamics of the Cardassian universe including the OC, 
MPC and EC models, for recent works on the Cardas- 
sian universe, see |19l - l22| . In the last section, we will give 
some discussions and conclusions. 



II. BRIEFLY REVIEW ON THE UNIFIED 
FIRST LAW 



Einstein equations can be written in a form called " uni- 
fied first law" based on the general definition of black 
hole dynamics on trapping horizon, which was proposed 
by Hayward and developed by Cai et al In 
the rest of this section, we will give a brief review on this 
unified first law in the 3 -I- 1-dimensional spherical sym- 
metric spacetime, in which the metric could be locally 
written in the double-null form as 



(1) 



where dfi^ is the line element of the 2-sphere with unit ra- 
dius, r and / are functions of {C^ , (T)- So, each symmet- 
ric sphere has two preferred normal directions, namely 
the null directions d/d£,^, which will be assumed future- 
pointing in the following. And also, we will assume the 
spacetime is time-orientable. The expansions of the ra- 
dial null geodesic congruence are defined by 



9± = 2r-^d±r, 



(2) 



where d± denotes the coordinates derivative along 
The expansion measures whether the light rays normal 
to the sphere are diverging (9± > 0) or converging 
(0^ < 0), namely, whether the sphere is increasing or 
decreasing in the null directions. Note that, although 
the value of 6± will change with geometries, its sign 
will not, and the only invariants of the metric and its 
first derivative are functions of r and e-^0+0_, or equiv- 
alently g°'^dardbr = —^6^9+9-, which has an important 
physical and geometrical meaning: a sphere is said to 
be trapped (untrapped) if 9+9^ > O(0+0_ < 0), and if 
9+9- = 0, it is a marginal sphere. 

Considering non-stationary black holes, Hayward has 
proposed that the future outer trapping horizons defined 
as the the closure of a hypersurface foliated by future 
or past, outer or innner marginal sphere is taken as the 
definition of black holes, since the horizon possess various 
properties which are often intuitively ascribed to black 
hole including confinement of observers and analogues 
of the zeroth, first and second law of thermodynamics. 
However, in the case of FRW universe, one should take 
the future inner trapping horizon defined by 



0, 



9-<0, d-9+>0, 



(3) 



as a system on which the thermodynamics will be estab- 
lished, since the surface gravity is negative on the cosmo- 
logical horizon. 



In a spherical symmetric spacetime, one can obtain the 
total energy inside the sphere with radius r by calculating 
the Minsner- Sharp energy given by 



r 

2G 



which is a pure geometric quantity and has much bet- 
ter properties than the other definitions of energy when 
one consider the case of non-stationary spacetime. The 
relation between the Minsner- Sharp energy and others 
could be found in ref. [3| • There are also two invariants 
constructed from the energy- momentum tensor T'^'^: 



W 
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T" ab 
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which are called the work density, and 5* called the energy 
fiux vector (also called the energy-supply vector), whose 
components are 



(6) 



Here and in the following, a, b denotes the two dimen- 
sion space normal to the sphere. With the help of the 
definition of Minsner-Sharp energy and the above two 
quantities, one can write the (0, 0) component of Ein- 
stein equations as a "unified first law" : 



dE :^A^ + WdV . 



(7) 



rich information, e.g. by projecting the unified first law 
along the trapping horizon, we can obtains the first law 
of black hole thermodynamics,which has the form Q 



{dE, z) 



where k defined by 
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(dA, z) + iy^dV, z) , 



1 



(8) 



(9) 



is the surface gravity of the trapping horizion. Here z — 
-f 2:~9_ is a vector tangent to the trapping horizon, 
and it should be noticed that by definition of the horizon 
d+r — 0, one has 

z''da{d+r)^z+d+d+r + z-d-d+r^Q (10) 

on the trapping horizon, then 

d+d+r 



z 



d-d: 



(11) 



Also note that by taking the Einstein equations = 
— 47rrr++, see ref. 3 and the definition of the surface 
gravity in Eq. ([9|), one can easily finds 



{A^,z) 



SttG 



{dA, z) , 



(12) 
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which is the Clausius relation in the version of black hole 
thermodynamics, see the first term on the right side of 
Eq. ([8]). The left side of the above equation is nothing 
but the heat flow 5Q, and the right side has the form 
TdS, if one identities the temperature T — n/iir and the 
entropy S — A/AG. So, in Einstein theory, the "uni- 
fied first law" also implies the Clausius relation, and this 
relation is also hold in the Gauss-Bonet and Lovelock 
gravity theories by treating the higher derivative terms 
as an effective energy-momentum tensor. But, in the 
scalar-tensor theory, this relation is no longer hold due 
to some non-equilibrium thermodynamical properties. In 
the next section, we will study the thermodynamics of the 
Cardassian universe, and calculate the corresponding en- 
tropy. 



III. THERMODYNAMICS OF THE 
CARDASSIAN UNIVERSE 

The spacetime of the Cardassian universe is described 
by the FRW metric, which could be written in the form 
of 



ds^ = habdx°'dx + fdfl^ 



(13) 



where = t, — r and f — a(t)r, which is the radius 
of the sphere while a{t) is the scale factor. Defining 



dS,^ = ^[dtT 



dr , (14) 



^/2 V VI - kr^ 



where k is the spacial curvature, the metric could be 
rewritten as a double-null form 



ds^ = -2d£.+dr + f-^dfl^ , 



(15) 



then we get the trapping horizon vt by solving the equa- 
tion d+f\f=f^ = as 



-1/2 



(16) 



where g is some function of the energy density of mat- 
ter and we have defined the effect energy density pe — 
g{pm) — Prn- Using the continuity equations, we can ob- 
tain the effective pressure corresponds to the effective 
energy density 



Pe = {Prn + Pm)g'iPm) ~ g{Pm) 



Pn 



(19) 



where the prime denotes the derivative with respect to 
Pm- Then, we get the associated work density We and 
energy-supply vector as 



We 



2g - Pm +Pm~ {Pm + Pm)g' 



(20) 



*e = -{g' -l)(p,n+Prn){-HrTdt + adr). (21) 
Therefore, we obtain 

SQ. . (.4*., .:^(!^).t(<!^) m .) 



2ttG V 9' 



where we have used the relation 



4G.g' 



(22) 



H--!^ = -^^~A7:Gg'iprn+Pm)- (23) 

Here, we have also identified T = k/2tt. And, for the 
heat flow of pure matter, we also have 

T 

{dA,z)-{A^e,z) = ^{dA,z), (24) 
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So, when gm = Pm, the above equation reduces to the 
Clausius relation in the unmodified Friedmann model. 
In the following, we will focus on some concrete Car- 
dassian models, namely, the original Cardassian model 
(OC), the modified polytropic Cardassian model (MFC) 
and the exponential model (EC). In these models, the 
function g{pm) takes different forms, which will reduce 
to Pm in the early universe and differ from the unmod- 
ified Friedmann universe at redshift z < 0(1), during 
which it will gives rise to accelerated expansion. 



which has the same form of the apparent horizon. Thus, 
one can get the surface gravity k = — (1 — ej/fr, where 
we have defined 



2HrT 



(17) 



One can also check that 9_fT < indicating the trapping 
horizon is future. By using Eq. (|10l) and after a direct 
calculation, one can get z~ = e/(l — e) when z"*" = 1 is 
chosen. Then, in the (t, r) coordinates, the project vector 
is given by z = 9t — (1 — 2e)Hrdr- 

In the Cardassian universe, the Friedmann equation is 
modified as 



k 8ttG . , SttG, , 

— = ^-5(Pm) = ^-(Pm +Pe) , (18) 



OC model 



In this model, the function g{pm) is given by 



aiPrn) = Pn 



Pm[l + foiPrn)] , (25) 



where pc is a character energy density in the Cardassian 
universe and the parameter n is assumed to satisfy n < 
2/3 to give rise to a acceleration of the universe. Here we 
have defined the function foiPm) = {pm/ Pc)^~^ ■ Thus, 
by using Eq. ([24| . the heat flow of pure matter is given 
by 



5Qn 



T 



4G(1 + nfo) 



{dA, z) = T{dSm, z) 



(26) 
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where the entropy is obtained by 



4G(1 + nfo) 



(27) 



which reduces to the usual relation dS = dA/AG in the 
limit of /o — 0. However, when fo is large enough, the 
entropy becomes a constant that independent of the sur- 
face area, which means there is no heat flow of pure mat- 
ter on the trapping horizon, but this time dSe — dA/AG. 
From Friedmann equation (HH]), we also have the follow- 
ing constraint equation 



fo-' (1 + fo) = 



2GAp, 

Therefore, we get the entropy by integrating Eq. (P7)) 
3 



(28) 



8G2pc(n-l) 



fo 



4G V 



2, 



-'{i + for^dfo 

-1 n-2 



n — 1 ' 71 — 1 ' 



-fo 



.(29) 



up to some integration constant. Here 2-F1 is the hyperge- 
ometric function and Eq. (|28)) gives the relation between 
fo and A. To illustrate the relation between the entropy 
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FIG. 1: Top: The entropy of the original Cardassian universe 
as the function of the surface area with parameter n — 0.1. 
Bottom: The area of the trapping horizon with respect to 
function fo with the same parameters. 



and the surface area, we plot an example with parameter 
n = 0.1 in Fig.[Tl in which it shows that when fo is small. 



the entropy satisfies the usual area law S = A/ AG, while 
fo is large, it becomes a constant as we see in Eq. (P7| . 
Actually, from Eq. (P^. one can obtain that 



n - 2\ ^ f2n - 1 



(30) 



where F denotes Gamma functions. 

B. MPC model 



The modified polytropic Cardassian model can be ob- 
tained by introducing an additional parameter q > 
into the original Cardassian model, in which the func- 
tion g{prn) is given by 



1+ ( ^ 



Pm[l + /q(Pm)]' 



(31) 

where fq — {pm/pc)'^^"' ^\ and when q = 1, it reduces 
to the original model. Thus, by using Eq. (|24p . the heat 
flow of pure matter is given by 



SQni = 



T 



4G(l + n/,)(l + /,)^ 
where the entropy is obtained by 

dA 



-{dA,z)^T{dS„ 



dS„ 



4G(l + n/,)(l + /,)^ 



-1 ' 



(32) 



(33) 



which also reduces to the usual relation dS — dA/AG in 
the limit of fq — > 0, and the entropy becomes a constant 
that independent of the surface area when fq is large 
enough. From Friedmann equation (jlSp . we also have 
the following constraint equation 



2GAp 

Therefore, we get the entropy by integrating Eq. ([33 

^ '' ^'^'\i + f,r^dfq 

1 ' 



(34) 



8G^p,qin~l) ' 



A 



4G^ + ^^ 



2^1 



-1 



9' ?("-!)' q{n-l) 



(35) 



up to some integration constant. Here Eq. (|34|) gives 
the relation between fq and A. To illustrate the rela- 
tion between the entropy and the surface area, we plot 
an example with parameter n = 0.1 in Fig. [2l in which 
it shows that when / is small, the entropy satisfies the 
usual area law S = A/ AG, while / is large, it becomes a 
constant as we see in Eq. ([55)) . Also, from Eq. (|55|) . one 
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FIG. 2: Top: The entropy of the original Cardassian universe 
as the function of the surface area with parameter n = 0.1 
and q = 0.5(solid), 0.7(dotted), l.O(dashed). Bottom: The 
area of the trapping horizon with respect to function fq with 
the same parameters. 



can obtain that 



r 1 



9( 



q{n-l) 



r 



(i) 



(36) 

so, when q = 1, it reduces to Eq. (15(1)) . Also, the top figure 
in Fig. [5] indicates that for a given n, all the curves cross 

the point 2^^^ , where / > satisfies / (1 + /) = 
1. 



C. EC model 

In this model, the function g{pm) is given by 



g{pm) = Pm exp 



Pm 
Pc 



Pme 



(37) 



where feipm) = (Pm/Pc) Again, by using Eq. 
the heat flow of pure matter is given by 

^ {dA,z) ^T{dS,n,z) , (38) 



4Ge^(l-n/,) 
where the entropy is obtained by 

dA 



dS„ 



4Ge/Kl-n/e) 



(39) 



which also reduces to the usual relation dS = dA/AG in 
the limit of fq — > 0, and the entropy becomes a constant 
that independent of the surface area when fq is large 
enough. From Friedmann equation (|18|) . we also have 
the following constraint equation 



2GApc ' 

Thus, the area has a maximum value 

3 



2Gp, 



■{ne)- 



(40) 



(41) 



when /e = 1/n. Therefore, we get the entropy by inte- 
grating Eq. ((39)) 



SG^PcU 
3 



-^^'f'e ^dfe 



8G2pcn2i/" 



r(l/7i)-r(l/n,2/e) 



, (42) 



up to some integration constant. Here Eq. (j^U)) gives 
the relation between /e and A. To illustrate the relation 




FIG. 3: Left: The entropy of the original Cardassian universe 
as the function of the surface area with parameter n = 0.6 
and the black point corresponds to the maximum point of the 
area. 



between the entropy and the surface area, we plot an 
example with parameter n = 0.6 in Fig. [3J in which it 
shows that when / is small, the entropy satisfies the usual 
area law S = A/AG, and it reaches its maximum value 



6 



at A = Ajnax when /e = 1/n. Actually, from Eq. (|42]) . 
one can obtain that 

8GV^'" ^ 4G' ^^^^ 
where we have used the relation pp)) . 

IV. DISCUSSION AND CONCULSION 

In this paper, we have studied the thermodynamics of 
the Cardassian universe and calculated its corresponding 
entropy, in particular, for the OC, MPC and EC Car- 
dassian models and the thermodynamic law is regarded 
as the origin of Cardassian model. We find that the en- 
tropy obeys ordinary area law on the trapping horizon 
when the area is small, and it becomes a constant when 
area is going to be large in the OC and MPC model, 
while it has a maximum value in the EC model. 

As we known that, the holographic principle states 
that all the information in the bulk should be encoded in 
the boundary of the bulk, so it seems that the Cardas- 
sian universe could only contain finite information. This 
may lead to a question that does the information in our 



universe will be infinite or not? Of course, if the ordi- 
nary area law S = A/ {AG) is always valid, the entropy 
will blow up when A goes to infinite value, in contrast 
with the case of that in the Cardassian models, S finally 
becomes a constant or gets its maximum value, when A 
is large. We will make a further study on this interesting 
topic [23*1. 

It should be noticed that, the cosmological constant 
will not change the area law S = A /{AG), because only 
the derivative g' emergences in Eq. p4l) or p2|). but 
in general, different dynamic dark energy models will 
change the law differently. So, the method we used in 
this paper provides a new way to distinguash different 
kinds of dark energy models. 
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